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Abstract. We propose an extension of a result by Repetowicz et al. 
1 about Wick's theorem and its applications: we first show that Wick's 
theorem can be extended to the uniform distribution on the sphere and 
then to the whole class of elliptical distributions. Then, as a special case, 
we detail this result for distributions that are Gaussian scale mixtures. 
Finally, we show that these results allow to recover easily a theorem by 
Folland 5 about integration of polynomials over the sphere. 



1. Wick's theorem 

In quantum field theory [2], the determination of the partition function 
of some systems involves the computation of integrals of the form 

where P (X) is a monomial in the components of I £ R n , that is, in a 
probabilistic approach, mixed moments of a Gaussian random vector with 
covariance matrix £. 

Wick's theorem provides a simple formula for these moments: let us de- 
note by X a Gaussian vector in W 1 with zero mean and covariance matrix 
S; a linear form C% on M. n is identified with the vector of its coefficients 

U (X) = ijx 

A pairing a on the set 12m = {1> • • • > 2m} is a partition of l 2 m into m disjoint 
pairs; the set Ii2 m of pairings of I 2m is composed of 2 2 ™m\ elements. For 
example, in the case m = 2, there are 3 different pairings on I\ = {1, 2, 3, 4} , 
namely a l = {(1, 2) , (3, 4)} , a 2 = {(1, 3) , (2, 4)} and a 3 = {(1, 4) , (2, 3)}. 
We denote by hm/v the set of indices i such that a = {(i, a (*))} • 
Wick's theorem states as follows [2]. 

Theorem. Let m be an integer then 

2m 

e\\c 1 (x)= n ^« 

If 2m is replaced by 2m + 1, then the expectation above equals zero. 

As an example, with m = 2 and Ci (X) = X; L , we obtain 

E [A 1 A 2 A 3 A 4 ] = E [XiX 2 ] E [X 3 A 4 ]+£ [X1A3] E [A 2 A 4 ]+£ [X1X4] E [X 2 X 3 ] . 

1 
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2. Spherical and Elliptical distributions 

In this paper, we extend Wick's theorem to the class of spherical random 
vectors: such a vector X G W 1 is characterized [31 Ch.2] by the fact that 
its characteristic function ip (Z) = Eexp [iZ l X\ is invariant by orthogonal 
transformation; in other words, 

^{z) = <j ) {z t z) 

for some function <j> : R + — > C. In this case, a stochastic representation of 
X writes 

X = rU 

where U is a random vector uniformly distributed on the unit sphere S n = 
{X £ R n ; ||X|| = 1} and r 6 M + is a positive random variable independent 
of U. If moreover vector X admits a density fx , then it is necessarily of the 
form 

fx (X) = g {X l X) 

for some function g : R — > R + . 

An extension of the spherical property is the elliptical property [31 Ch.2]: 
a random vector Y £ R n is elliptical with (n x n) characteristic matrix X if 

(1) Y = T} /2 X 

where X is spherical. We assume without loss of generality that all vectors 
have zero mean. 

3. Wick's theorem for the uniform distribution on the sphere 

Wick's theorem can be extended to the case of a vector uniformly dis- 
tributed on the sphere S n as follows. 

Theorem. // U is uniformly distributed on the sphere S n then 

g n^>= 2mr( l?U s n 

Proof. A stochastic representation for vector U is 

v v 



\N\\ 

where iV is a Gaussian vector with unit covariance matrix. By the polar 
factorization property, random variable r = ||jV|| is independent of U; we 
deduce from the linearity of Ci that 

2m 2m 2m 

E\[C i {N) = E\[rC i (U) = Er 2m E [] d (U) 

i=l i=l i=l 

so that 

i=i 
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Since r = \/N f N is chi-distributed with n degrees of freedom, elementary- 
algebra yields 

, 2m _ 2"T(m + f) 

r(f) 

□ 



Er l 



From the linearity of forms £,% and by stochastic representation (pQ), we 
deduce a version of Wick's theorem for any elliptical vector as follows. 

Theorem. If X £ R n is an elliptical vector with stochastic representation 

(2) X = aT, 1/2 U 

and if random variable a has moments up to order 2m, then 

e u^ x )=w^e n ^« 

4. Wick's theorem for Gaussian scale mixtures 

A random vector X is a Gaussian scale mixture if its probability distri- 
bution writes 



r+oo 

fx(X) = / g(x;aZ)dH(a) 
Jo 



where g (x; E) is the Gaussian distribution with (n x n) covariance matrix 
S on W 1 and H is a cumulative distribution function. A stochastic repre- 
sentation of such a vector is 

(3) X = ^ 

where iV is a Gaussian vector with covariance matrix E and a is a scalar 
random variable independent of N with cumulative distribution function H. 
Cauchy distributions, alpha-stable and Student-t distributions belong to the 
family of Gaussian scale mixture distributions. 

If X is a Gaussian scale mixture, then it is also elliptical; however, the 
contrary is not true, as shown [4] by the Pearson II (or Student-r) family of 
distributions 



fx(X) 



r (f=f + i) (_ x t K- 1 x\i-<> 



rfe)|7rpJf|Va 



p 



with p = n+ 2y^ > and notation {x) + = max (x, 0) . 

Using stochastic representation ©, Wick's theorem extends to Gaussian 
scale mixtures as follows. 
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Theorem. If X is a Gaussian scale mixture as in with Ea~ m < oo, 
then 

2m 

E\{C l {X) = E[a~ m } IT ^» 

i=l crgn m i£l2 m /a 

We note that Repetowitz et al. [lj deal only with the case of Student-t 
distributions for which random variable a is Gamma distributed; the pro- 
posed theorem extends thus their result to the class of all Gaussian scale 



5. Application: a result by Folland 

In this section, we show that the results above allow to recover easily the 
following theorem by Folland [5j about integration of a monomial over the 
sphere S n . 

Theorem. // 

n 

then 



i=l 



where 



f JO if some oii is odd 

/ P (X) da = < 2T{fr)...T{p n ) j 

a °i + 1 



2 

and da is the surface measure on the sphere S n . 

Proof. Defining «o = 0, we consider the function T : [l,X^=i a *l — ¥ t 1 ' 71 ] 
such that 

T (i) = j if a\ + ■ ■ ■ + ctj-i < i < a\ + • • • + ctj. 

We choose forms £{ such that li = S^u) where 5k is the k-th column vector 
of the n x n identity matrix. Since the scalar product l\l a u\ equals zero 
when T (i) ^ T [a (i)] , we consider only the Y\^=i ^f a ^ Pairings for which 

r(i) = 6 T{i) 6 T(i) 



T (i) = T [a (i)] and = SL r \8 T u\ = 1. Thus, by Wick's theorem, 



2"T(m + f)^ 2 ^(|)! 2-r(m+§) (2^^, 



Since da is the unnormalized surface measure and as the surface of the 
sphere is S (<S n ) = |^y, 

r. n 

/ P(X)da = S(S n ) xET] Ci{U) 
which yields the result. □ 
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